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CLAIMS: 

1 . A method of deriving blood perfusion indices for a region of interest (ROI) of a 
subject, the method comprising the steps of: 

5 administering a contrast agent to the subject during a dynamic imaging scan: 

converting signal intensity data from raw images of the scan into contrast agent 
concmtration data; 

deriving parameters from the contrast agent concentration data using at least one 
transport function that accotmts for delay and dispersion of the contrast agent; and 
10 calculating the blood perfusion indices from the derived parameters. 

2. A method according to claim 1 wherein the transport function represents a 
probability distribution frmction of transit times of the contrast agent through the 
subject 

15 

3. A method according to claim 2 further comprising the step of using a first model 
to represent an arterial transport function ha(t) through a vessel leading to the ROI. 

4. A method according to claim 3 further comprising using a second model to 
20 represent a tissue transport function hs(t) through the ROI. 

5. A method according to claim 4 further comprising the step of selecting an 
arterial input function AIFa(t) in the vessel leading to the ROI by searching pixels taken 
of the contrast ^ent concentration data. 

25 

6. A method according to claim 5 further comprising the step of measuring the 
contrast agent concentration C(t) remaining in the ROI. 

7. A method according to claim 6 further comprising the step of representing ha(t) 
30 using a gamma*-variate function (GVF) in the first model such tiiat: 
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where =CTj*'^r(l+ai), T (a ) s x*'*e"'dbc is the Gamma function, ti is 

the time taken for the contrast agent to move firom the initial measurement of AIFa(t) to 
a vessel at the entry to the ROI, ai and ai are related to the mean transit time and 
dispersion of ha(t). 

8. A method according to claim 7 further comprising the step of estimatiug ha(t) 
after deriving values for parameters ti and oi and setting =0 using the equation: 



0 (t<t,) 



10 9. A method according to claim 8 further comprising the step of determining an 
estimate for the arterial input function AIFt(t) of the vessel at the entry to flie ROI using 
the equation: 

where 0 is the convolution operator. 

15 

10. A method according to claim 9 further comprising the step of determining an 
estimate of blood flow Ft and an estimate of the tissue JKF Rc(t) jfrom the deconvolution 
of: 

C(0 = (Ft/kidAIF^t) 0Re(t) 

20 

where kH==(l-Ha)/(l-Ht) is a correction constant taking into account different values of 
arterial hematocrit Ha and tissue hematocrit Ht since the contrast agent remaius in the 
extracellular fraction of blood (plasma). 

25 11. A method according to claim 10 further comprising the step of determining an 
estimate for the tissue transport function he(t) from the estimated Re(t) using the 
equation: 

he(t)^-j^Re(t) 
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10 



12. A method according to claim 11 further comprising the step of determining a 
rise time and a mean transit time of hc(t) in order to determine parameters and <52 by 
assuming t2=0, where t2, aa and <5% are parameters related to the mean transit time and 
dispersion of he(t). 

13. A method according to claim 11 further comprising the step of determining a 
peak height and a mean transit time of he(t) in order to det^mine parameters oz and 1% 
by assuming a2=0, where t2, a2and 02 are parameters relating to mean transit time 
and dispersion of he(t). 

14. A method according to claim 12 or claim 13 further comprising the step of 
representing a simulated transport function hs(t) using a GVF in the second model such 
that: 



15 



where = a^^'^Ti} + ^2 ) , t2, 02 and az are parameters related to the mean transit time 
and dispersion of hs(t) through the ROI. 

15. A method according to claim 14 further comprising the step of estimating hs(t) 
20 using the derived values for parameters a2 and 02 by setting t2=0 using the equation: 



25 



16. A method according to claim 14 further comprising the step of estimating hs(t) 
using the derived values for parameters a2 and tiby setting using the equation: 



1 



,-(l-<j)/«r. 



^,(0 = 
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17. A method according to claim 15 or claim 16 fiirther comprising the step of 
determining a simulated tissue IRF RsCt) using the eqxiation: 
J?,(0-l-JjA,(r)c/tr 

5 18. A method according to claim 17 further comprising the step of det^mining a 
simulated contrast agent concentration Cs(t) using the equation: 
Q(t)^ (Ft/kid AIF^t)0R^(t) 

19. A method according to claim 18 further comprising the step of fitting the 
10 simulated Cs(t) to C(t) using a least squares method according to: 

20. A method according to claim 19 further comprising the step of optimising the 
parameters Ft, ti,ai, ai, 02, Qq, and t2by minimizing S iteratively. 

15 

21. A method according to claim 20 further comprising the step of reducing the 
number of adjustable parameters by fixing ai=0 and t2=0, or fixing ai=0 and a2=0 
leading to five adjustable parameters. 

20 22. A method according to claim 20 or claim 21 comprising the step of further 
reducing the number of adjustable parameters by fixing a relative dispersion, 
Pi=ai/(ai+ti), of ha(t) resulting in ai dcpmd&at on t^^ and therefore leading to four 
adjustable parameters. 

25 23. A method according to claim 22 further comprising the step of calculating 
quantitative blood perfusion indices firom the optimized parameters of Fu ti, oi, ai, ozy 
a2 and t2. 

24. A method according to claim 23 wherdn the perfusion indices include any one 
30 or more of blood flow, blood volume, mean transit time, arterial delay time, arterial 
dispersion time or relative arterial dispersion, tissue dispersion time or relative tissue 
dispersion. 
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25. A method according to claim 24 fbrfher comprising tiie step of repeating eadi 
previous step, apart from the step of selecting the AIF, on a pixel-by-pixel basis to 
produce quantitative maps of the perfusion indices for further analysis and presentation. 

5 26. A method according to any one of claims 1 to 25 wherein the ROI is a tissue. 

27. A method according to any one of claims 1 to 25 wherein the ROI is a pixel or a 
plurality of pixels in a tissue. 

10 28. A method according to any one of claims 1 to 27 wherein the scan is any one of 
CT,MRIorNM. 

29. A method according to any one of claims 3 to 24 wherein the vessel is an artery. 

15 30. A method according to claim 29 further comprising determining a venous input 
function VIFa(t) from a draining vein to estimate an AIFa(t) where a selected artery has 
partial voluming, the vein being larger than the artery. 

31. A method according to claim 30 further comprising the step of determining the 
20 profile ofVIFa(t) from the draining vein. 

32. A method according to claim 3 1 further comprising the step of scaling AIFa(t) to 
have the same first-pass bolus peak area as the VIFa(t) to minimize partial voluming 
effect from the AIFa(t). 



33. A method according to claim 32 wherein the first-pass bolus peak areas of the 
AIFa(t) and VIFa(t) profiles are obtained by fitting the profiles to ganoma-variate 
function (GVF) profiles respectively to remove contrast recirculation effects. 

30 34. A method according to any one of claims 17 to 33 further comprising the step of 
detenxuning a simulated tissue IRF Rs(t) in the case that the contrast agent does not 
always remain in the vascular system, such as in a tumour in the subject in order to 
determine blood perfusion indices and permeability indices using: 



25 




35 
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-(f-/2)/<rj 



where 



0 



E is the extraction fraction of the tracer in the blood stream that leaks out of the vessel 
into tissue, and the tracer clearance rate constant k=E*Ft A^'e is a rate constant at which 
the leaked contrast agent diffuses back into the blood stream and leaves the tissue, Ve is 
5 volume fraction of the extravascular and extracellular space (EES). 

35. A method according to claim 33 wherein a permeability sur&ce area product PS 
is detOTnined hy PS = -F^ hi(l - E) . 

10 36. Computer program means for deriving blood perfusion indices for a region of 
interest (ROT) of a subject by directing a processor to carry out any of the method steps 
according to any one of claims 1 to 35 apart from the step of administering a contrast 
agent to the subject during a dynamic imaging scan. 

15 37. Computer program means according to claim 36 further directing the processor 
to retrieve raw image data from the dynamic imaging scan of the subject after a contrast 
agent is administered to the subject. 

38. A system of deriving blood perfusion indices for a region of interest (ROI) of a 
20 subject, the system comprising: 

scanning means for providing a dynamic image scan of the subject during which 
a contrast agent is administered to the subject; 

processor means linked to the scanning means for retrieving raw image data 
from the scan; 
25 the processor means further. 

converting signal intensity data included in the retrieved raw image data into 
contrast agent concentration data; 

deriving parameters from the contrast agent concentration data using at least one 
transport function that accounts for delay and dispersion of the contrast agent; and 
30 calculating the blood perfusion indices from the derived parameters. 
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39. A system according to claim 38 wherein the transport fimction represents a 
probability distribution function of transit times of the contrast agent through the 
subject 

S 40. A system according to claim 39 wherein a first model is used to repTesent an 
arterial transport function ha(t) through a vessel leading to the ROI. 

41. A system according to claim 40 wherein a second model is used to represent a 
tissue transport function hs(t) through the ROI. 

10 

42. A system according to claim 41 wherein the processor means selects an arterial 
input function AIFa(t) in the vessel leading to the ROI by searching pixels taken of the 
contrast agent concentration data. 

15 43. A system according to claim 42 wherein the processor means measures the 
contrast agent concentration C(t) remaining in the ROI. 



